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ABSTRACT 

The Cygnus Loop was observed from the northeast to the southwest with 
XMM-Newton. We divided the observed region into two parts, the north path 
and the south path, and studied the X-ray spectra along two paths. The spectra 
can be well fitted either by a one-component non-equilibrium ionization (NEI) 
model or by a two-component NEI model. The rim regions can be well fitted by 
a one-component model with relatively low kTe whose metal abundances are sub- 
solar (0.1-0.2). The major part of the paths requires a two-component model. 
Due to projection effects, we concluded that the low kT^ (~0.2keV) component 
surrounds the high kT^ (~0.6keV) component, with the latter having relatively 
high metal abundances (~5 times solar). Since the Cygnus Loop is thought to 
originate in a cavity explosion, the low kT^, component originates from the cavity 
wall while the high kTg component originates from the ejecta. 

The flux of the cavity wall component shows a large variation along our path. 
We found it to be very thin in the south-west region, suggesting a blowout along 
our line of sight. The metal distribution inside the ejecta shows non-uniformity, 
depending on the element. O, Ne and Mg are relatively more abundant in the 
outer region while Si, S and Fe are concentrated in the inner region, with all 
metals showing strong asymmetry. This observational evidence implies an asym- 
metric explosion of the progenitor star. The abundance of the ejecta also indicates 
the progenitor star to be about 15 Mq. 

Subject headings: ISM: abundances — ISM: individual (Cygnus Loop) — super- 
nova remnants — X-rays: ISM 
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Introduction 



A supernova remnant (SNR) reflects the abundance of tlie progenitor star wfien tfie 
remnant is young and tliat of tlie interstellar matter (ISM) when it becomes old. In this 
way, we can study the evolution of the ejecta and the ISM. The Cygnus Loop is a proto- 
typical middle-aged shell-like SNR. The angular diameter is about 2°. 4 and it is very close 
to us (540 pc; Blair et al. 2005), implying a diameter of ~23pc. The estimated age is 
about 10000 year s, less than half that based on the previous distance estimate of 770 pc 
(iMinkowski IQSsh . 



Since the Cygnus Loop is an evolved SNR, the bright shell mainly consists of a shock- 
heated surrounding material. Its supernova (SN) explosi on is generally considered to have 
occurred in a preexisting cavity flMcCray fc Snow 19791 ). Levenson et al. (1997) found 
that the Cygnus Loop was a result of a cavity explosion that was created by a star no 
later than BO. It is almost circular in shape with a break-out in the south where the hot 
plasma extends out of the circular shape. Miyata et al. (1994) observe d the northeast (NE ) 
shell of the Loop with ASCA and revealed the metal deficiency there (IMiyata et al. 19941 ). 
Since Dopita et al. (1977) reported the metal deficiency of the ISM around the Cygnus 
Loop, they concluded that the plasma in the NE-shell is dominated by the ISM. Due to 
the constraints of the detector efficiency, they a ssumed that the relative abundances of C, 
N and O are equal to those of the solar valu e (lAnders fc Grevesse 19891 ). More recently, 
Miyata et al. (2007) used the Suzaku satellite ( Mitsuda et al. 20071 ) to observe one pointing 
position in the NE rim. They detected emission lines from C and N and determined the 
relative abundances (IMiyata et al. 20071 ). They concluded that the relative abundances of C, 
N and O are consistent with those of the solar values whereas the absolute abundances show 



depletion from the solar values (lAnders &: Grevesse 19891 ). Katsuda et al. (2007) observed 
four pointings in the NE rim and detected a region where the relative abundances of C and 
N are a few times higher than that of O. 

Hatsukade & Tsunemi (1990) detec ted a hot plasma inside the Cygnus Loop that is 
not expected in the simple Sedov model (IHatsukade fc Tsunemi 19901 ). They reported that 
the hot plasma was confined inside the Loop. Miyata e t al. (1998) detect ed strong emis- 
sion lines from Si, S and Fe-L from inside the Loop flMiyata et al. 19981 ). They found 
that the metal abundanc e is at least several times higher than that of the solar value 
( lAnders &: Grevesse 19891 ). indicating that a few tens of higher than that of the shell re- 
gion. They concluded that the metal rich plasma was a fossil of the SN explosion. The 
abundance ratio of Si, S and Fe indicated the progenitor star mass to be 25 M©. Miyata 
& Tsunemi (1999) measured the radial profile inside the Loop and found a discontinuity 
around 0.9 Rg where Rg is the shock radius. They measured the metallicity inside the hot 
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cavity and estimated the progenitor mass to be 15 Mq. Levenson et al. (1998) estimated the 
size of the cavity and the progenitor mass to be 15 Mq. Therefore, the progenitor star of the 
Cygnus Loop is a massive star in which the triple-a reaction should have dominated rather 
than the CNO cycle. If the surrounding material of the Cygnus Loop is contaminated by 
the stellar activity o f the progenitor star , it may explain the C abundance inferred for this 



region with Suzaku (IKatsuda et al. 20071 ). 



In order to study the plasma condition inside the Cygnus Loop, we observed it from the 
NE rim to the south-west (SW) rim with the XMM-Newton satellite. We report here the 
result covering a full diameter by seven pointings. 



2. Observations 

We performed seven pointing observations of the Cygnus Loop so that we could cover 
the full diameter from the NE rim to the SW rim (from Pos-1 to Pos-7) during the AO-1 
phase. We concentrate on the data obtained with the EPIC MOS and PN cameras. All the 
data were taken by using medium filters and the prime full window mode. Fortunately, all the 
data other than Pos-4 suffered very little from background flares. Obs IDs, the observation 
date, the nominal point, and the effective exposure times after rejecting the high-background 
periods are summarized in table [TJ 

All the raw data were processed with version 6.5.0 of the XMM Science Analysis 
Software (XMMSAS). We selected X-ray events corresponding to patterns 0-12 and for 
MOS and PN, respectively. We furt her cleaned t he data by removing all the events in 



bad columns listed in the literature (IKirsch 20061 ). After filtering the data, they were 



vignetting-corrected using the XMMSAS task evigweight. For the background subtrac- 
tion, we employed the data set accumulated from blank sky observations prepared by 



(iRead &: Ponman 20031 ). After adjusting its normalization to the source data by using the 



energ y range between 5 keV and 12key , where the emission is free from the contamina- 
tion (IFujita et al. 2004 : ISato et al. 20051 ). we subtracted the background data set from the 



source. 



-4- 



3. Spatially Resolved Spectral Analysis 
3.1. Band image 

Figure [T] displays an exposure-corrected ROSAT HRI image of the entire Cygnus Loop 
(black and white) overlaid with the XMM-Newton color images of the merged MOSl/2, PN 
data from all the XMM-Newton observations. In this figure, we allocated color codes as red 
(0.3-0.52 keV), green (0.52-1.07keV) and blue (1.07-3 keV). We see that the outer regions 
are reddish rather than bluish while the central region is in bluish. 

The NE rim is the brightest in our field of view (FOV), showing a bright filament at 45° 
to the radial direction correspondin g to NGC6992. The SW ri m is also bright in our FOV 



where there is a V-shape structure (lAschenbach fc Leahy 19991 ). In the center of the Loop 



an X-ray bright filament runs through Pos-4 and Pos-5 forming a circular structure. In the 
ROSAT image, we can see it and find that it forms a large circle within the Cygnus Loop. 
In this way, there are many fine bright filaments in intensity. However, we find that there is 
a clear intensity variation along our scan path: dim in the center and bright in the rim. 

Figure [2] shows spectra for seven pointings; each is the sum of the entire FOV. The NE 
rim (Pos-1) and the SW rim (Pos-7) show strong emission lines below IkeV including O, 
Fe-L and Ne, while the center (Pos-4) shows strong emission lines from Si and S. We can see 
that the equivalent width of Si and S emission lines are bigger in the center and gradually 
decrease toward the rim. We show the comparison of spectra between Pos-1 and Pos-4 in 
figure O Prominent emission lines are 0-Hea, 0-Lya, Fe-L complex, Ne-Hea, Mg-Hea, 
Si-Hea, and S-Hea. We see that the emission line shapes for O are quite similar to each 
other while there is a big difference at higher energy band. Since t he spectrum in the NE 
rim can be well represented by a single temperature plasma model (jMiyata et al. 1994 ). we 
need an extra component in the center. 



3.2. Radial profile 

Although there are many fine structures, no matter how finely we divide our FOV, each 
region would contain different plasma conditions due to the integration of the emission along 
the line of sight. Therefore, we concentrate on large scale structure along the scan path. 
First of all, we divided our FOV into two parts along the diameter: the north path and the 
south path. Then we divided them into many small annular sectors whose center is located 
at (20''51™34^7, 31°00'00"), i.e., the center of the nominal position of Pos-4. There are 141 
and 172 annular sectors for the north path and the south path, respectively. These small 
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annular sectors, shown in figure [T], are divided such that each has at least 60,000 photons 
(~20,000 for MOSl/2 and ~40,000 for PN) to equalize the statistics. We extracted the 
spectrum from each sector using the data set accumulated from blank sky observations as 
sky background. We have confirmed that the emission above 3 keV is statistically zero. In 
this way, we obtained 313 spectra. These sectors can be identified by the angular distance, 
"R", from the center (east is negative and west is positive as shown in figure [1]). 

The width of the sector depends on R. The sector widths range from 3'. 8 to 0'.2 in the 
north path and from 3'.0 to 0'.2 in the south path. The widest sectors are in Pos-4 due to 
its short exposure because of the background flare. The narrowest sectors are in the NE rim 
where the surface brightness is the highest. 



3.3. Single temperature NEI model 



We fitted the spectrum for each sector wit h an absorbed non-equilibrium ionization 
(NEI) model with a single kTp, using the Wabs (IMorrison fc McCamnaon 19831) and VNEI 
model (NEI version 2.0) JSorkowski et al. 200lh in XSPEC v 12.3.1 JArnaud 1996hl. We 
fixed the column density, iVn to be 4.0x 10^°cm"^ (e.g., (Ilnoue et al. 19801 ). (IKahn et al. 1980l )) 
Free parameters were kTe, the ionization time scale, r (a product of the electron density 
and the elapsed time after the shock heating); the emission measure (hereafter EM; EM = 
/ neUndi, where nn and Ue are the number densities of hydrogens and electrons, di is the 
plasma depth); and abundances of C, N, O, Ne, Mg, Si, S, Fe, and Ni. We set abundances 
of C and N equal to that of O, that of Ni equal to Fe, and other elements fixed to the 
solar values (lAnders fc Grevesse 19891 ). In the fitting process, we set 20 as the minimum 
counts in each spectral bin to perform the test. We determined the value of the minimum 
counts such that it did not affect the fitting results. Figure H] shows the distribution of the 
reduced-x^ in black as a function of R along both the north path and the south path. We 
found that values of reduce d-y^ for all the sectors are between 1.0 and 2.0. If we took into 
account a systematic error (INevalainen et al. 20031 ; iKirsch 20041 ) of 5%, the reduced-x^ was 
around 1.5 or less. 

In general, values of the reduced-x^ are a little higher in the central part of the Cygnus 
Loop. Miyata et al. (1994) observed the NE rim with ASCA and found that the spectra were 
well represented with a one-temperature VNEI m odel with a temper ature gradient towards 
the inside. The Suzaku observation in the NE rim (IMiyata et al. 20071 ) reveals that the X-ray 
spectrum can be represented by a two-temperature model: one component is 0.2-0.35 keV 
and the other is 0.09-0.15 keV. In our fitting, the value of kTe obtained is 0.2-0.25 keV. 
Therefore, we detected a hot component that Suzaku detected. There may be an additional 
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component with low temperature that seems difficult to detect with XMM-Newton due to 
the relatively lower sensitivity below 0.5 keV compared to Suzaku. 



The ASCA observation (IMiyata et al. 19941 ) also shows that the metal abundance in 



the NE rim is deficient. The authors concluded that the plasma in the NI3-rim consists of 
the interstell ar matter (ISM) r ather than the ejecta. This is confirmed with the Suzaku 



observation (iMiyata et al. 20071 ) that indicates the abundances of C, N and O to be ~0.1, 
0.05 and 0.1 solar , respectively. We also obtained the metal deficiency in the NE rim data; 
the best-fit results are given in figure [5] (left) and in table O Leahy (2004) measured the 
X-ray spectrum of the southwest regi on of the Cy gnus Loop and reported that the oxygen 



abundance there is about 0.22 solar (ILeahy 20041 ). Therefore, the X-ray measurements of 



the Cygnus Loop show that the metal abundances are depleted. 

Cartledge et al. (2004) measured the interstellar oxygen along 36 sight lines and con- 
firmed the homogeneity of the 0/H ratio within 800 pc of the Sun. We found that they 
measured it in the direction about 5° away from t he Cygnus Loop. The oxy gen abundance 



they measured is about 0.4 times the solar value (lAnders fc Grevesse 19891 ). Wilms et al. 
(2000) employed 0.6 of the total interstellar abundances for the gas-phase ISM oxygen abun- 
dance, and suggest that this depletion may be due to grains. Although the ISM near the 
Cygnus Loop may be depleted, the abundances are still much higher than what we obtained 
at the rim of the Cygnus Loop. It is difficult to explain such a low abundance of oxygen 
in material originating from the ISM. Therefore, the origin of the low metal abundance is 
open to the question. Since the Cygnus Loop is thought to have exploded in a pre-existing 
cavity, we can say that the cavity material shows low metal abundance. The abundance 
difference between our data and those from Suzaku may be due to the difference in the 
detection efficiency at low energy. Taking into account the projection effect, the plasma of 
the rim regions consists only of the cavity material while that of the inner regions consists 
both of the cavity material and of an extra component filling the interior of the Loop. 



3.4. Two-temperature NEI model 

To further constrain the plasma conditions, we applied a two-component NEI model 
with different temperatures. In this model, we added an extra component to the single 
temperature model. The extra component is also an absorbed VNEI model with kT^, r, 
and EM as free parameters. The metal abundances of the extra component are fixed to 
those determined at the NE rim so that the extra component represents the cavity mate- 
rial. Figure H] shows reduced-x^ values in red along the path. Applying the F-test with 
a significance level of 99% to determine whether or not an extra component is needed, we 
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found that most of the spectra required a two-component model, particularly in the central 
part of the Loop. Sectors that do not require two-component model are mainly clustered 
in R<-65', +25' <R< +40', and +60' <R. Therefore, we considered that the outer sectors 
(|R|> 70') can be safely represented by a one-component model while other sectors can be 
represented by a two-component model. In this way, we performed the analysis by applying 
a two-component VNEI model with different temperatures. We assumed that the low tem- 
perature component comes from the surrounding region of the Cygnus Loop and the high 
temperature component occupies the interior of the Loop. 

We found that the values of reduced-x^ are 1.0-1.8 even employing a two-component 
model. This is partly due to the systematic errors. Looking at the image in detail, there 
are fine structures within the sector. Furthermore, the spectrum from each sector is an 
integration along the line of sight. Since we only employ two VNEI plasma models, the 
values of reduced-x^ are mainly due to the simplicity of the plasma model employed here. 
Therefore, we think that the plasma parameters obtained will represent typical values in 
each sector. 

Figure (right) and table [3] shows an example result that comes from the sector at 
R=+10'. Fixed parameters in the low fcT^ component come f rom the fitting result at the 



NE rim obtained by Suzaku observations (lUchida et al. 20061 ). Metal abundances for the 



high kTe component show higher values by an order of magnitude than those of the low 
fcTg component, surely confirming that the high fcTg component is dominated by fossil ejecta. 

Figure [6] shows temperatures as a function of position. The low kT^ component is in 
the temperature range of 0.12-0.34 keV while the high kTe component is above 0.35 keV. 
There is a clear temperature difference where a two-component model is required rather 
than a single temperature model. The low kT^ component represents the cavity material 
surrounding the Cygnus Loop while the high fcTg component represents the fossil ejecta 
inside the Loop. Figure [7] shows the fluxes for the two components as a function of position. 
The low kTe component shows clear rim brightening. The east part is stronger than the west 
part, showing asymmetry of the Loop. On the other hand, the high fcTg component has a 
relatively flat radial dependence. From the center to the SW, we see that the flux of the 
high kTe component is stronger than that of the low kT^ component. 



3.4.1. Distribution of the cavity material 



As shown in flgureO the low kTe component shows relatively constant temperature with 
radius. The distribution of the flux shown in flgure [7] shows peaks in the rim and relatively 
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low values inside the Loop. There are some differences between the north path and the south 
path. The biggest one is a clear difference in peak position in the NE rim that is due to the 
bright filament at 45° to the radial direction, as seen in figure [H However, these two paths 
show a globally similar behavior in flux. Therefore, we can see that they are quite similar to 
each other from a large scale point of view. 

We notice that there are many aspects showing asymmetry and non-uniformity. The 
NE half is stronger in intensity than the SW half. The flux in the inner part of the Loop 
shows relatively small values in the west half, particularly at +25' <R< +40'. The NE half 
is brighter by a factor of 5-10 than the SW half. Furthermore, the SW half shows stronger 
intensity variation than the NE half. This suggests that the thickness of the cavity shell is 
far from uniform. The cavity shell in the SW half is much thinner than that in the NE half. 
Since we assumed the metal abundances of the low kT^ component equal to those of the NE 
rim, we can calculate the EM. Furthermore, we assumed the ambient density to be 0.7 cm~^ 
based on the observation of the NE rim, and we estimate the mass of the low kT^ component 
to be 130 Mq. However, we should note that there are is evidence that the SN explosion 
which produced the Cygnus Loop occurred within a preexisting cavity (e.g., Hester et al. 
1994; Levenson et al. 1998; Levenson et al. 1999). The model predicts that the original 
cavity density, Uc, is related to the wall density Ug by ric = 5. Assuming that no equals the 
ambient density, n^, we estimate ric to be 0.14 cm~^. Then, we calculate the total mass in 
the preexisting cavity to be ~25 M©. 

3.4-2. Ejecta distribution 

The flux distribution from the ejecta along the path is shown with filled circles in figure[71 
It has a relatively fiat structure with two troughs around R=— 35' and R=+50'. Since we 
left the metal abundances as free parameters, we obtained distributions of EM of various 
metals (C=N=0, Fe=Ni, Ne, Mg, Si and S) in the ejecta. These are shown in figure[8], where 
black crosses trace the north path and red crosses trace along the south path. If we assume 
uniform plasma conditions along the line of sight, the EM represents the mass of the metal. 
Most elements show similar structure between the north path and the south path while there 
is a big discrepancy in Fe/Ni distribution at —10' < R < +30'. In this region, the south 
path shows two times more abundant Fe/Ni than the north path. A similar discrepancy is 
seen in O (-30' < R < -10') and in Ne (at -10' < R < +10'). Therefore, the distribution 
of metal abundance shows a north-south asymmetry along the path. 

The distributions of O and Ne show a central bump and increase at the outer sectors. 
However, those of Mg, Si, S and Fe only show a central bump. The increase of O and Ne 
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at the outer sectors indicates that the outer parts of the ejecta mainly consists of O and Ne 
and they may be well mixed. Similarly heavy elements, Mg, Si, S and Fe/Ni forming central 
bumps may show that they are well mixed. Therefore, a significant convection has occurred 
in the central bumps while an "onion-skin" structure remains in the outer sectors. 

4. Discussion 

The Cygnus Loop appears to be almost circular with a blowout in the south. The 
ROSAT image indicates no clear shell in this blowout region. Levenson et al. (1997) revealed 
that there is a thin shell left at the edge of the blowout region. Therefore, there is a small 
amount of cavity material in this region that surrounds the ejecta. This also indicates the 
non-uniformity of the cavity wall. If the cavity wall is thin, the ejecta can produce a blowout 
structure. 

Looking at the component of the cavity material along our path shown in figure [71 the 
flux is very weak at -|-15' <R< +40'. This indicates that the cavity wall is very thin in this 
region. When we calculate the flux ratio between the ejecta plasma and the cavity material, 
we find that the ratio becomes high (larger than 4) at +15' <R< +35' in the north path 
and +30' <R< +35' in the south path. Therefore, we guess that the thin shell region is 
larger in the north path than in the south path. This also shows the asymmetry between 
the north and the south as well as that between the east and the west. If the thin shell 
region corresponds to a blowout similar to that in the south blowout, this region must have 
a blowout structure along the line of sight either in the near side or far side or both. This 
structure roughly corresponds to Pos-5 and will extend further in the northwest direction. 
Looking at the ROSAT image in figure [H we see a circular region with low intensity. It is 
centered at (20^49'"11'', 31°05'20") with a radius of 30'. We guess that this circular region 
corresponds to a possible blowout in the direction of our line of sight. CCD observation just 
north of our path will answer this hypothesis. 

We obtained EMs of O, Ne, Mg, Si, S, and Fe for the ejecta along the north path and the 
south path. Multiplying the EMs by the area of each sector, we obtained emission integral 
(hereafter EI, EI= J neUudV, dV is the X-ray-emitting volume) along the path. Since we 
only observed the limited area of the Cygnus Loop from the NE rim to the SW rim, we 
have to estimate the Els for the entire remnant in order to obtain the relative abundances 
as well as the total mass of the ejecta. Therefore, we divide our observation region into four 
parts: left- north part, right-north part, left-south part, and right-south part. We assume 
that each part represents the average Els of the corresponding quadrant of the Loop. In 
this way, we can calculate the total Els for O, Ne, Mg, Si, S, and Fe that are described in 
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table m The south quadrant, corresponding to the right-south path, contains the largest 
mass fraction of 31%, while the other quadrants contain 23% each. Then, we calculate the 
relative abundances of Ne, Mg, Si, S, and Fe to O in the entire ejecta. Since we cannot 
measure the abundance of light elements like He, it is quite difficult to estimate the absolute 
abundances. However, the relative abundance to O is robust. 

Since the Cygnus Loop is believed to be a result from a core-collapse SN, we compared 
our data with core-collapse SN models. There are many theoretical results from various au- 
thors (e.g., Woosley & Weaver 1995; Thielemann et al. 1996; Rauscher et al . 2002; Tominaga 



et al. 2007). We also employed a SN Type la model (llwamoto et al. 19991 ) for comparison. 
We calculated the relative abundance for various elements to O and compared them with 
models. Figure M shows comparisons between the model c alculations and our result s where 



we picked up Woo sley's model with one sola r abundance (I Anders fc Grevesse 19891 ) for the 



core-collapse case (jWoosley fc Weaver 19951 ). The Type la model yields more Si, S and Fe 
than our results, but less Ne. Models with massive stars produce better fits to our results 
than the Type la model. Among them, we found that the model with 15 Mq showed good 
fits to our results. They fit within a factor of two with an exception of Fe. We also noticed 
that the model with one solar abundance looked better fit than that with depleted abun- 
dance. Therefore, we can conclude that the Cygnus Loop originated from an approximately 
15 Mq star with one solar abundance. 

Assuming that the ejecta density is uniform along the line of sight, we estimate the 
total mass of the fossil ejecta to be 21 Mq. In this calculation, we assumed that the electron 
density is equal to that of hydrogen and that the plasma filling factor is unity, although 
the fossil ejecta might be deficient in hydrogen. If it is the case, the total mass of the 
fossil ejecta reduces to ~12, Mq whereas the relative abundances are not affected. The most 
suitable nucleosynthetic model predicts that the total mass ejected is about 6 Mq without H. 
Therefore, there might be a significant amount of contamination from the swept-up matter 
into the high-fcTe component, which we consider the ejecta. Otherwise, the assumption that 
the density of the ejecta is uniform might be incorrect since rim-brightening for the EMs of 
O, Ne, Mg, and Fe is clearly seen in Fig. [HI Non- uniformity reduces the filling factor and 
also the mass of the high-fcTe component. 

There is observational eviden ce of the asymm etry of supernova exp losions both for 



massive stars (ILeonard et al. 20061 ) and for Type la (IMotohara et al. 20061 ). We found that 
the ejecta plasma shows asymmetric structure between NE half and SW half. Ne and Fe are 
evenly divided while two thirds of O and Mg are in the NE half. On the contrary, two third of 
Si and S are in the SW half. We calculated the ejecta mass for each quadrant and found that 
the south quadrant contains the largest ejecta mass. Similar asymmetries are seen in other 
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SNR, such as Puppis A, which shows asymmetric structure with 0-rich, fast-moving knots 
(Winkler & Kirscher 1985; Winkler et al. 1988) . The central compact object in Puppis A 
is on the opposite side of the SNR from the 0-rich, fast-moving knots (jPetre et al. 19961 ). 
If the asymmetry of the ejecta in the Cygnus Loop is similar to that of Puppis A, we may 
expect a compact object to be in the north direction. 



5. Conclusion 

We have observed the Cygnus Loop along the diameter from the NE rim to the SW rim 
employing XMM Newton. The FOV is divided into two paths: the north path and the south 
path. Then it is divided into many small annuli so that each annulus contains a similar 
number of photons to preserve statistics. 

The spectra from the rim regions can be expressed by a one-fcTg component model while 
those in the inner region require a two-kT^ component model. The low kT^ plasma shows 
relatively low metal abundance and covers the entire FOV. It forms a shell that originates 
from the preexisting cavity. The high kTe plasma shows high metal abundance and occupies a 
large part of the FOV. The origins of these two components are different: the high kTe plasma 
with the high metal abundance must come from the ejecta while low kT^ plasma with low 
metal abundance must come from the cavity material. We find that the thickness of the 
shell is very thin in the south west part where, we guess, the ejecta plasma is blow out in 
the direction of our line of sight. 

We estimate the mass of the metals. Based on the relative metal abundance, we find that 
the Cygnus Loop originated from a 15 Mq star. The distribution of the ejecta is asymmetric, 
suggesting an asymmetric explosion. 
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of Education, Culture, Sports, Science and Technology (16002004). This study is also carried 
out as part of the 21st Century COE Program, ^Towards a new basic science: depth and 
synthesis'. S. K. is supported by JSPS Research Fellowship for Young Scientists. 



- 12 - 

REFERENCES 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Arnaud, K. A. 1996, ASP Conf. Ser., 101, 17 

Ashcehnbach, B. & Leahy, D. A. 1999, A&A, 341, 602 

Burrows, Adam; Hayes, John, 1996, Physical Review Letters, 76, 352 

Blair, W. P., Sankrit, R., & Raymond, J. C., 2005, AJ, 129, 2268 

Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. 2001, ApJ, 548, 820 

Cartledge, S. I. B., Lauroesch, J. T., Meyer, D. M., & Sofia, U. J. 2004, ApJ, 613, 1037 

Dopita, M. A., Mathewson, D. S., and Ford, V. L. 1977, ApJ, 214, 179 

Fujita, Y., Sarazin, C. L., Reiprich, T. H., Andernach, H., Ehle, M., Murgia, M., Rudnich, 
L., & Slee, O. B. 2004, ApJ, 616, 157 

Ghavamian, P., Raymond, John., Smith, R. C., & Hartigan, P. 2001, ApJ, 547, L995 

Hatsukade, I., & Tsunemi, H. 1990, ApJ, 362, 566 

Hester, J. J., Raymond, J. C., & Blair, W. P. 1994, ApJ, 420, 721 

Inoue, H., Koyama, K., Matsuoka, M., Ohashi, T., Tanaka, Y., Tsunemi, H. 1980 ApJ, 238 
886 

Iwamoto, Brachwitz, F., Nomoto, K., Kishimoto, N., Umeda, H., Hix, W. R., & Thielemann, 
F. K., 1999, ApJS, 125, 439 

Kahn, S. M., Charles, P. A., Bowyer, S., Blissett, R. J. 1980, ApJ, 242, L19 

Katsuda, S., Tsunemi, H., Uchida, H., Miyata, E., Nemes, N., Miller, E. D., and Hughes, J. 
P. 2007, PASJ in press 

Kirsch, M. 2004, EPIC status of calibration and data analysis, XMM-SOC-CAL-TN-0018, 
issue 2.3 

Kirsch, M. 2006, XMM-EPIC status of cahbration and data analysis, XMM-SOC-CAL-TN- 
0018, issue 2.5 

Leahy D. A. 2004, MNRAS, 351, 385 



Leonard, D. C, Filippenko, A. V., Ganeshalingam, M., Serduke, F. J. D., Li, W., Swift, B. 
J., Gal- Yam, A., Foley, R. J., Fox, D. B., Park, S., Hoffman, J. L., and Wong, D. S., 
2006, Nature, 440, 505. 

Levenson, N. A., et al., 1997, ApJ, 484, 304. 

Levenson, N. A., Graham, J. R., Keller, L. D., Richter, M. J. 1998, ApJS, 118, 541 
Levenson, N. A., Graham, J. R., & Snowden, S. L. 1999, ApJ, 526, 874. 438, LI 15 
McCray, R. & Snow, T. R, Jr. 1979, ARA&A, 17, 213. 
Minkowski, R., Rev. Mod. Phys. 1958, 30, 1048. 
Mitsuda, K. et al., 2007, PASJ, 59S, 1. 

Miyata, E., Tsunemi, H., Pisarki, R., and Kissel, S. E. 1994, PASJ, 46, LlOl 

Miyata, E., Tsunemi, H., Kohmura, T., Suzuki, S., and Kumagai, S. 1998, PASJ, 50, 257 

Miyata, E., & Tsunemi, H. 1999, ApJ, 525, 305 

Miyata, E., Katsuda, S., Tsunemi, H., Hughes, J. P., Kokubun, M., & Porter, F. S. 2007, 
PASJ, 59S, 163. 

Morrison, R. & McCammon, D. 1983, ApJ, 270, 119 

Motohara, K., Macda, Keiichi, G., Christopher L., Nomoto, K., Tanaka, M., Tominaga, N., 
Ohkubo, T., Mazzah, P. A., Fesen, R. A., Hoflich, P., Wheeler, J. G. 2006, ApJ. 652 
LlOl 

Nevalainen, J., Lieu, R., Bonamente, M., & Lumb, D. 2003, ApJ, 584, L716 
Petre, R., Becker, C. M., & Winkler, P. F. 1996, ApJ, 456, L43 
Rauscher, T., Heger, A., Hoffman, R. D., Woosley, S. E. 2002, ApJ, 576, 323 
Read, A., M., & Ponman, T. J. 2003, A&A, 409, 395 

Sato, K., Furusho, T., Yamasaki, Y., Ishida M., Matsushida, K., and Ghashi, T. 2005, PASJ, 
57, 743 

Thielemann, F.-K., Nomoto, K., & Hashimoto, M., 1996, ApJ, 460, 408 
Tominaga, N., Umeda, H., Nomoto, K. in prep. 



-14- 



Uchida, H. Katsuda, S., Miyata, E., Tsunemi, H., Hughes, J. P., Kokubun, M., & Porter, F. 
S. 2006, Suzaku Conference in Kyoto 

Wilms, J., Allen, A., and McCray, R. 2000, ApJ, 542, 914 

Winkler, P. F. & Kirshner, R. P. 1985, ApJ, 299, 981 

Winkler, P. F., Tuttle, J. H., Kirshner, R. P., & Irwin, M., J. 1988, in lAU CoUoq. 101: 
Supernova Remnants and the Interstellar Medium, ed. R. S. Roger & T. L. Landecker, 
65-+ 

Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 15 - 



Table 1: Summary of the seven observations. 



Obs. ID 


Camera 


Obs. Date 


Coordinate (RA, DEC) 


Effective Exposure 


0082540101 


MOSl 


2002-11-25 


20'*55'"23^6, 31°46'17".0 


14.1 ksec 


(Pos-1) 


M0S2 






14.1 kscc 




PN 






5.6kscc 


0082540201 


MOSl 


2002-12-03 


20''54"07^4, 31°30'51.4".0 


14.4 ksec 


(Pos-2) 


M0S2 






14.4 ksec 




PN 






11.7 ksec 


0082540301 


MOSl 


2002-12-05 


20'^52"^51M, 31°15'25".7 


11.6 ksec 


(Pos-3) 


M0S2 






11.6 ksec 




PN 






9.1 ksec 


0082540401 


MOSl 


2002-12-07 


2o/»5im34. 7^ 31°00'00".0 


4.9 ksec 


(Pos-4) 


M0S2 






4.9 kscc 




PN 






3.4 ksec 


0082540501 


MOSl 


2002-12-09 


20''50'"18".4, 30°44'34".3 


12.6 ksec 


(Pos-5) 


M0S2 






12.6 ksec 




PN 






10.0 ksec 


0082540601 


MOSl 


2002-12-11 


20M9"'02^0, 30°28'16".l 


11.5 ksec 


(Pos-6) 


M0S2 






11.5 ksec 




PN 






5.9 ksec 


0082540701 


MOSl 


2002-12-13 


2o/»47m45s 8^ 30°13'42".9 


13.5 ksec 


(Pos-7) 


M0S2 






13.5 ksec 




PN 






7.5 ksec 
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Table 2: Spectral-fit ])aiaiiioters. 



Paraiiietor 


region— 7 1.25 


ArH[1020cm-2 


4 (fixed) 


A;Te[keV] 


0.23 ±0.01 


0(=C=N) 


0.068 ±0.002 


Ne 


0.17±0.01 


Mg 


0.14±0.03 


Si 


0.3±0.1 


S 


0.6±0.2 


Fe(=Ni) 


0.157±0.006 


log(T/cm-3sec) . 11.31 ±0.02 


EM^fxlO^^ cm-5] ll-Olas 


xVd-oi 


420/314 



Note. Other elements are fixed to solar values. 
The values of abundances are multiples of solar value. 
The errors are in the range Ax^ < 2.7 on one parameter. 
-■^EM denotes emission measure, J neUuM. 
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Table 3: Spectral-fit parameters. 



Parameter 


region+10 


A^H[102°cm-2 


4 (fixed) 



Low temperature component 



A;re[keV] 0.20 ±0.01 

C 0.27 (fixed) 

N 0.10 (fixed) 

O 0.11 (fixed) 

Ne 0.21 (fixed) 

Mg 0.17 (fixed) 

Si 0.34 (fixed) 

S 0.17 (fixed) 

Fe(=Ni) 0.20 (fixed) 

log(r/cm"''' sec) . 12 < 

EM^[xlO^^ cm-5] 1.34l^°;^o4 

High temperature component 

A;re[keV] 0.48 ±0.01 

0(=C=N) <0.01 

Ne 0.15l[i:[|? 

Mg 0.21±0.08 

Si 2.5±0.3 

S 5±1 

Fe(=Ni) 1.03±0.04 

log(r/cm-3 sec) . . 11.12±0.05 

EM^[xl0^9 cm-5] 0.094+°:°°^ 

xVd.o.f 531/377 



Note. Other elements are fixed to solar values. 
The values of abundances are multiples of solar value. 
The errors are in the range Ax^ < 2.7 on one parameter. 
^EM denotes emission measure, / Uenudi. 
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Table 4: Calculated emission integral of the Cygnus Loop ejecta. 



Element 10^'^ cm ^ 






7.4±0.5 


Ne 


1.5±0.2 


Mg 


0.34±0.1 


Si 


2.9±0.5 


S 


1.2±0.3 


Fe 


1.30±0.05 




Fig. 1. — Left: Exposure-corrected ROSAT HRI image of the entire Cygnus Loop (black 
and white) overlaid with the XMM- Newton color images of the merged MOSl/2, PN data 
from all the XMM-Newton observations. Right: Spectral extraction regions overlaid on the 
XMM-Newton three-color image shown in left figure. 
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Fig. 2. — MOSl spectra for the seven pointings; each is the sum of the entire FOV. 
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Comparison of spectra between Pos-1 (circles) and Pos-4 (triangles). They are 
in intensity at O-Hea. 
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Fig. 4. — Distribution of the reduced-x^ as a function of R along the north path (upper 
panel) and the south path (lower panel). The single-component model is shown in black and 
the two-component model is shown in red. 
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Fig. 5. — Left: An example spectrum that comes from the sector at R=— 74.25'. The best-fit 
curves are shown with sohd hues and the lower panels show the residuals. Right: Same as 
left but for the sector at R=+10'. Both the ejecta and cavity components are shown only 
for MOSl spectrum as dashed lines. 
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Fig. 6. — Temperature distributions of the two components as a function of position. Filled 
circles show the ejecta component, while crosses show the cavity component. Black show 
the north path and red shows the south path. Typical errors are ±5% for both components. 
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Fig. 7. — Flux distributions of the two components as a function of position. The marks in 
this figure are the same as those in figure El 
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Fig. 8.— Distributions of EM for various metals (0(=C=N), Ne, Mg, Si, S, and Fe(=Ni)) 
in the ejecta. Black indicates the north path and red indicates the south path. Data points 
showing only upper limits are excluded in these figures. 
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Fig. 9. — Number ratios of Ne, Mg, Si, S, and Fe relative to O of the high-/cTe component 
estimated for the entire Loop (black sohd hn e). Dash-dotted red lines represent CDDl and 
W7 model of Type-la (llwamoto et al. 19991 ) . Dotted green, blue, hght blue, and magenta 
lines represer it core- collapse models w hose progenitor masses are 13, 15, 20, and 25 M©, 
respectively (jWoosley fc Weaver 19951 ). 



